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ABSTRACT 

We present the star formation history and chemical evolution of the Sextans 
dSph dwarf galaxy as a function of galactocentric distance. We derive these 
from the VI photometry of stars in the 42' x 28' field using the SMART model 
developed by Yuk & Lee (2007, ApJ, 668, 876) and adopting a closed-box model 
for chemical evolution. For the adopted age of Sextans 15 Gyr, we find that 
>84% of the stars formed prior to 1 1 Gyr ago, significant star formation extends 
from 15 to 11 Gyr ago (~ 65% of the stars formed 13 to 15 Gyr ago while ~ 
25% formed 11 to 13 Gyr ago), detectable star formation continued to at least 8 
Gyr ago, the star formation history is more extended in the central regions than 
the outskirts, and the difference in star formation rates between the central and 
outer regions is most marked 11 to 13 Gyr ago. Whether blue straggler stars are 
interpreted as intermediate age main sequence stars affects conclusions regarding 
the star formation history for times 4 to 8 Gyr ago, but this is at most only a 
trace population. We find that the metallicity of the stars increased rapidly up to 
[Fe/H]=-1.6 in the central region and to [Fe/H]=-1.8 in the outer region within 
the first Gyr, and has varied slowly since then. The abundance ratios of several 
elements derived in this study are in good agreement with the observational data 
based on the high resolution spectroscopy in the literature. We conclude that the 
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primary driver for the radial gradient of the stellar population in this galaxy is 
the star formation history, which self-consistently drives the chemical enrichment 
history. 

Subject headings: galaxies: stellar content — galaxies: chemical abundance — galaxies:evolution 
- galaxies: individual(Sextans) 



1. Introduction 



We know of abou t a dozen dwarf spheroidal (dSph) satellite galaxies of the Milky Way 



Galaxy (IMateo I Il998l ; iBellazzini et al.ll2006l ; iTolstoy et al.l 120091 ) . Although these galaxies 



are among the faintest and least massive, they are, judging from the discovered population 
in the Local Group, the most common in the Universe. As such, they provide interesting 
constraints on models of galaxy formation and evolution in a cosmological context. The 
nearest of these, because of the increased sensitivity and angular resolution, provide the 



galaxies. In turn, such studies inform investigations of the first stars in the Universe ( 


Wvse 


2005; 


Ricotti & Gnedin 


2005: 


Kawata et al. 


2006; 


Fenner et al. 


2006: 


Bromm et al. 


2009) 



and the stellar halo of the Milky Way Galaxy 

For various reasons, the SFH of these galaxies was originally presumed to consist of a 
simple single burst of star formation at early times. In part, this was the natural extension of 
the known characteristics of previously identified halo populations (the old stellar halo and 
the old, single burst population of globular clusters). However, deep photometric studies 
found that the SFH of dSph's is much more complex and intriguing than expe cted (e.g, 
Lee et ali[l~993l : IMateo Ifl998l : bolpin II2OO2I : bolpin et al1l2005l : llblstov et al]l2009h and have 
ultimately shown that simple star formation histories are the exception rather than the rule. 



990). is one of the fainter known 



(120031 ) presented deep, wide field 



The Sextans dSph, discovered in 1990 (llrwin et al. 
galaxies (My > — 12 mag) in the Local Group. iLee et al. 
VI photometry of this galaxy based on the CFH12 K images, and d iscussed inferences drawn 
from the color-magnitude diagram. To summarize, ILee et al.l (120031 ) derived a schematic star 
formation history of the Sextans dSph and found the following: 1) most stars in Sextans 
formed at a similar time, and with similar low metallicity, as the metal-poor globular, M92, 
2) hints of continuing star formation in the inferred metallicity spread (A[Fe/H]~ 0.2 dex) 
among red giant branch stars, population gradients, strong red horizontal branch, and bright 
mai n-sequence stars sl i ghtly above the main sequence turnoff of the main population (see 
also llkuta fc Arimoto I (120021 )) and 3) a small fraction of all the stars, including anomalous 
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Cepheids, formed less than 1 Gyr ago. 



Somewhat in contrast, the large metallicity range o f the individual s tars ( —1.45 < [Fe/H]< 
—2.85), derived for five red giant stars in Sextans by IShetrone et al.l (120011 ) from high res- 
olution spectra obtained with the High Resolution Echelle Spectrometer (HIRES) on the 
Keck I telescope is inconsistent with the hypothesis of a dominant, short-lived star forma- 
tion episode. The existence of such a metallicity spread differs from what is observed in 
globular clusters and requires some self-enrichment, which in turn implies an extended SFH. 



To resolve the potential discrepancy among various studies, we reanalyze the iLee et al. 



(2003) photometr y using the SMART model, a new numerical algorithm developed by 
Yuk fc Lee I (120071 ). which quantitatively recovers the SFH and chemical evolution by synthe- 
sizing the color-magnitude diagram. We discuss the following as follows: in §2 the SMART 
model used to derive the SFH and chemical evolution, and the available data for Sextans, 
in §3 the results of the modeling, in §4 comparisons to other models and the implications of 
the results, and in §5 our summary of the main results. 



2. Model and Data 
2.1. Model 

The SMART model recovers the SFH and chemical evolution of a stellar system by 
synthesizing the color-magnitude diagram in comparison with observations. The principal 
difference between the numerical algorithm in SMART and other algorithms used to derive 
the SFH of resolved stellar populations is that SMART treats the in ternal chemical evo lution 



self-consistently. The detailed description of this model is given in lYuk fc Lee I (120071 ). and 
so we describe it briefly here. The core aim of this model is to simultaneously derive the SFH 
and chemical evolution of galaxies from observational color-magnitude diagrams (CMDs). 
The algorithm produces a set of simulated CMDs for varying star formation histories, and 
compares these to the observed CMD to fin d the best-fitting SF H. In this model we use 



the Padova database of theoretical isochrones flGirardi et al. Il2002l ). The Padova isochrones 



cover the evolution of stars from zero-age main sequence to carbon ignition for higher mass 
stars and to the thermally pulsating asymptotic giant branch phase for lower mass stars. 
We preform interpolation to make a fine grid using logarithmic variables according to the 
evolutionary equivalent point of evolutionary tracks. 

The best-fit is determined by minimizing the goodness-of-fit parameter x 2 - If the best- 
fit star formation rate at the k-th time bin is 4> s k ° l , then the asymmetric confidence limits of 
the fc-th time bin can be obtained as follows. We find the best-fit star formation rate with 
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the fixed <fi s k ° l + 5<fik, and calculate goodness-of-fit parameter + ^4>k)- If the difference 

between Xpi^t" 1 ) an d xl(4>T l + ^4>k) is lcr, then is the la confidence limit of the fc-th time 
bin. In the similar way, we calculate confidence limits of each time bin in plus and minus 
direction. 



2.2. Data 



We u s e the VI photometry of stars in a 28' x 42' field centered on Sextans given by 



Lee et al.l (120031 ). This photometric catalog comes from CCD images obtained using the 
CFH12K camera at the CFHT. The limiting magnitudes (50% co mpleteness ) are V = 24.4 
mag and I = 23.6 mag. To investigate spatial variations in the SFH. lLee et al.l ( 120031 ) radially 
divided the entire field into four different regions: Rl(< 7'.3), R2(7'.3 - 11'), R3(ll' - 15'), 
and R4(> 15'). The boundaries of the regions are set to maintain a constant number of 
mid-red giant branch (RGB) stars with 0.8 < (V — I) < 1.1 mag and 18 < V < 21 mag. 
The properties of these regions are summarized in Table 1. For reference, the boundary of 
region R3 (15 arcmin) is slightl y smaller than the core radius of the Sextans dSph 16'. 6 ± 1'.2 
(jlrwin fc Hatzidimitriou Ill995l ). 



In Figure [T] we pr esent the observe d stellar CMDs in each of the four regions of the 
Sextans dSph given by iLee et al.l (120031 ). Bright foreground stars were not used in model- 
ing. The five r ed gia nt branch stars for which chemical abundances have been measured 
Shetrone et al.l (120011 ) are marked with circles. Basic information for these giant stars is 
presented in Table 2. 



One outstanding feature of this galaxy is its large number of blue stragglers (ILee et al. 



20031 ). which we highlight with pentagons in Figure 1. There are 59, 64, 68, and 69 such 
stars in Rl, R2, R3, and R4, respectively. 



3. Results 

We run the SMART algorithm in each of the four regions independently to investigate 
whether there are radial SFH variations. The model parameters we adopt for our calculation 
are summarized in Table [3] and we comment on a few of our choices. Previous estimates 
for the distance to Sextans range from (m — M)o = 19.7 to 19.9, and the reddening values 
range from E(B — V) — 0.01 to 0.03 (see the summary in Lee et al. ( 20031 )). Comparing the 



observational CMD and the synthetic CMD, we determined which distance produces the best 
model fits, obtaining a distance modulus (m — M) = 19.67 and reddening E(B — V) = 0.03 
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(E(V —I) = 0.04). This value of the distance is slightly smaller than that derived in lLee et al. 
(120031 ). At this distance one arcmin corresponds to 25.0 pc. Our choice of Sa lpeter's initial 
m ass func t ion (M m) = Am x ~ l ) with an index of x = —1.35 ( Salpeterl 1955 ). is supported 
by IWyse I (120051 ) who found that the initial mass function in Ursa Minor, a dSph similar to 
Sextans, is co nsistent this function . We adopted a closed- b ox m odel for chemical evolution 
model. While iFenner et al.l (120061 ) and llkuta fc Arimoto I (120021 ) suggested in the study of 
chemical evoluti on of dSphs that the galactic winds d o no t make a major role in removing 
gas from dSphs. lLanfranchi fc Matteucci I (120041 . 120071 ) and lRecchi et al. I (120081 ) presented a 
different view: chemical properties of the dSphs can be explained well only by the chemical 
evolution model involving intense galactic winds and low star formation efficiency In this 
study we will show that the closed box model also reproduces well the observed abundance 
ratios in Sextans. We also assume that metals produced in a region stay in that region (i.e., 
very inefficient mixing). We used the VI photometry brighter than the 60% completeness 
limit for modeling. We also assume that the current stellar mass of the galaxy was the same 
as the gas mass in the beginning. 



Tolstoy et al.l (120031 ) measured the ages of red giants with known abundances in this 



galaxy to be 15 Gyr using stellar models. This age is l arger than the age of the Universe 
derived from current high-precision cosmological models (IKomatsu et al.ll2009l ) and suggests 
that there is an unresolved problem with the stellar models. Matching the stellar ages of 
old populations with the age of the Universe has been a long-stand ing problem, and coul d 
even be considered the first evidence for a cosmological constant (IBolte fc Hogan I Il995l ); 
however, the current level of discrepancy is within the range of various plausible model 
adjustments. We adopt an age of 15 Gyr so that the current stellar models match the oldest 
populations, but caution against an absolute interpretation of the derived ages and favor a 
relative interpretation. We used seven age bins for calculation: 0-4, 4-7, 7-9, 9—11, 11-13, 
13-14, and 14-15 Gyr. 



3.1. Star Formation History 

We derive the SFHs of each of the four regions for the two blue-straggler scenarios (Case 
A and Case B). We normalized the star formation rate with respect to the lifetime average 
star formation rate for the entire observed region of Sextans, 1.7 x 10 -4 M Q yr _1 . In Figure[2] 
we present the observed CMD for the entire observed area and the best-fit synthetic CMD 
for Case A. The two CMDs look very similar in general, showing that the synthetic CMD 
reproduces well the observed CMD. We describe the resulting SFHs for Case A and Case B 
below. 



- 6- 



3.1.1. Case A : Including Blue Straggler Stars 

In Table H] we list the best-fit, Case A SFH and in Figure 3 we display it. The star 
formation rates are weighted by the number density of each region, and the errors correspond 
to the 1 cr confidence level. From the global SFH shown in Figure 3, we conclude that star 
formation was extended, and not a single burst like in Galactic globular clusters. While most 
star formation occurred between 15 and 7 Gyr ago, detectable star formation continued until 
4 Gyr ago. To be more quantitative, we find that 57 % of the stars formed between 13 and 
15 Gyr ago, 27 % of stars formed between 11 and 13 Gyr ago, the remainder formed mostly 
until 7 Gyr ago, and a tracer population, which consists of the blue stragglers, formed until 
4 Gyr ago. Thus, a majority (84 %) of the stars formed within the first 4 Gyr. 

The star formation rate decreased more rapidly in the outer regions. In the central 
region (Rl), 52 %, 32 %, and 10 % of stars formed 13 - 15, 11 - 13, and 9 - 11 Gyr ago, 
respectively, while in the outer region (R4), 72 %, 15 %, and 8 % of stars formed during the 
same periods. Thus, the star formation rate in the central region decreased by a factor of 
five from 14 Gyr to 10 Gyr ago, while that in the outer region decreased by a factor of nine. 
The most distinct period occurred 11-13 Gyr ago, when the star formation rate in the inner 
region was a factor of two larger than in the outer region. 

In Figure H] we display the radial behavior of the star formation rate within each age 
bin. There is an evident radial gradient in the star formation rate at all times when the star 
formation is significantly different than zero. Prior to 11 Gyr ago, the star formation rate 
decreases approximately linearly with increasing radial distance, and the gradient is steepest 
during the period 11-13 Gyr ago. Therefore, it is this population of stars, formed between 
11 and 13 Gyr ago, that is primarily responsible for the observed population gradient. 

3.1.2. Case B : Excluding Blue Straggler Stars 

The Case B SFH of each region is listed in Table [5] and is displayed in Figure [51 While 
the pattern of star formation is similar to the Case A results, one notable difference arises 
due to the exclusion of blue straggler stars in Case B: the Case B, SFH does not contain the 
tracer population that forms between 4-7 Gyr. 

Again, we find that the star formation decreases more rapidly in the outer regions. In 
the central region, Rl, 56 %, 34 %, and 7% of stars formed 13 — 15, 11 — 13, and 9 — 11 
Gyr ago, respectively, while in the outer region, R4, 81 %, 11 %, and 6 % of stars formed 
over the same periods. In total, we find that 65 % of the stars formed between 13 and 15 
Gyr ago, 25 % of stars were formed between 11 and 13 Gyr ago, and the remainder formed 
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until 7 Gyr ago. Again, we find that the vast majority (90%) of stars formed prior to 11 Gyr 
ago. This quantitative measurement varies only slightly between Case A (84%) and Case B 
(90%). 

In Figure [6] we display the radial variation of the star formation rate derived for Case B. 
As in Case A, the radial gradient of the star formation rate is prominent in age bins where 
we measure a non-zero star formation rate. The star formation rate decreases approximately 
linearly with increasing radial distance from the beginning, and the radial gradient is steepest 
11-13 Gyr ago. The shape of the radial distribution appears almost identical to that for Case 
A. 



3.2. Chemical Evolution: Metal Enrichment History 



The results for chemical evolution are very similar between Cases A and B so that we 
describe only the results for Case B hereafter. In Figure [7] we present the metal enrichment 
histories of each of the four regions in the Sextans dSph for Case B. We plot the mass- 
weighted mean metallicity of the stars versus time. The metallicity increased rapidly within 
~ 0.1 Gyr, and more slowly thereafter. 

Radially, the metallicity of the stars increased to [Fe/H]=-1.6 in the central region and 
to [Fe/H]=-1.8 in the outer region within the first Gyr. The early metal enrichment (within 
2 Gyr from the beginning) of Rl coincides with that of R2, because the star formation rate 
in Rl is nearly the same as that in R2. The metallicity difference between Rl and R4 was 
A[Fe/H] « 0.15 until 14 Gyr ago, and « 0.33 7 Gyr ago. 

In Figure [8] we display the radial metallicity distributions for four different ages (14.98, 
14.9, 14.0, and 7 Gyr). The radial gradient is zero within R = 10' and is -0.017 dex/arcmin 
at R > 10' within the first Gyr. The slope of the radial gradient becomes constant over the 
entire region (~ —0.03 dex/arcmin) 7 Gyr ago. 

Previous estima tes of the mean metallicity of stars in the Sextans dSph, as summarized 
by iLee et al.l (120031 ). were derived from the photometry of a large number of red giants, 
low-resolution spectroscopy, or high-resolution spectroscopy of a small number of red gi- 



ants. Interestingly, th e distribution of these estimates is rough 



at [F e/H] 1.6 dex (IMateo et all Il99ll : iDa Costa et al. 



1998 ), and the other at [Fe/H] 2.05 dex (ISuntzeff et al. 



1991 



y bimodal, with one peak 



Mateo et al 



199 



5; 



Mateo 



19931 ; iGeisler k, Sarajedini 



19961 : 



Shetrone et al.ll200ll ). The bimodality is present among both the photo metric an d spec- 



troscopic estimates. No reason for this bimodal distribution is known. ILee et al.l (120031 ) 



concluded on the basis of averaging that the mean metallicity of the giants in the Sextans 
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dSphs is ([Fe/H]) = —2.1 ± 0.2 with a dispersion of 0.2 dex. The metallicity estimates for 
five red giants based on high-resolut ion spectra show a l arge range from [Fe/H]=-2.85 to 
-1.45, with a mean of —2.07 ± 0.21 (IShetrone et al.ll200ll ). However, this is based on only 
five stars so we need to increase the sample size. 

In Table |6] we present the mean metallicities from our model. Our derived values of 
([Fe/H]) are independent of whether we include or exclude blue straggler stars. For our 
entire Sextans surve y region, we find that ([Fe/H ]) = —1- 7 6 dex , which is 0.3 dex larger than 



the values given by IShetrone et al.l (120011 ) and iLee et al.l (120031 ). However, this difference is 
insignificant (< 2a) because of the large dispersion among the observed [Fe/H] for individual 
stars and the small number of stars with spectroscopic measurements. 



3.3. Chemical Evolution: Chemical Abundance Ratios 

The evolution of the chemical abundance of a galaxy depends on its SFH because dif- 
ferent element s are produced d uring the evolution of stars with different masses over a range 



of timescales (IMatteucci II2008I ). A nice summa ry of the main featur es of the chemical evolu- 



tionary history of different elements is given by lTolstoy et al.l (120031 ). Most elements used in 
the study of chemical evolution of stellar populations belong to four groups: 1) light elements 
including O, Na, Mg, and Al, 2) a-elements (or even-Z element) including O, Mg, Si, Ca, and 
Ti, 3) Iron-peak elements including V, Cr, Mn, Co, Ni, Cu, and Zn, and 4) heavy elements 
including Y, Ba, Ce, Sm, and Eu. Light elements are a tracer of deep-mixing abundance 
patterns in RGB stars, and they show significantly different patterns between globular clus- 
ter giants and field stars, a-elements are mostly produced in massive stars and ejected in 
SNe II. Iron-peak elements are produced via the explosive nucleosynthesis sequence during 
the nuclear fusion (in both SN la and SNe II). Cu and Mn are considered to be produced 
by SNe la, in which case the ratio [Cu/a] or [Mn/a] is a good indicator of the ratio between 
SNe la and SNe II. Heavy elements provide a measure of the ratio of s-process (occurring in 
AGB stars) to r-process (occurring in SNe II) elements. 

A particularly informative parameter in the study of chemical evolution is the ratio 
[O/Fe]. Oxygen is produced primarily in high-mass stars of very short lifetimes and is 
ejected by SNe II, while iron is produced in both SNe II and SNe la. In Figure [9] we show 
the variation of [O/Fe] over time derived from Case B model. The value of [O/Fe] decreases 
very rapidly in the beginning, then more slowly after reaching [0/Fe]~ 0.08, and then even 
more slowly after 5 Gyr. [O/Fe] is similar to the solar value at about 11 Gyr ago. Radially, 
[O/Fe] varies little, except for 12 to 9 Gyr ago when it was slightly higher in the inner 
region than in the outer region. Unfortunately, it is difficult to observationally determine 
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[O/Fe] for stars, and the value of [0/H] is typically derived using emission line regions in a 
galaxy. Therefore, there are no available stellar [O/Fe] measurements in Sextans that can 
be compared to our predictions. 

Another informative parameter is the ratio [a/Fe]. This ratio, like [O/Fe], depends on 
the relative contribution of SNe II and SN la. Unlike [O/Fe] it can be measured from stellar 
spectra. The timescale for changes in [a/Fe] depends not only on the SFH, but also on other 
factors such as the initial mass function, the SNe la timescale, and the timescales for mix ing 
the SNe la and SNe II products back into the interstellar medium ( Matteucci 2003al lbh . 
We use a definition of the average [a/Fe] as an average of Mg , Ca , and Ti abundance s 
{[a/Fe] = [|[(Mg + Ca + Ti)/Fe]) following Ishetrone et al.l (boOlh and lTolstov et all (b()03h . 
In Figure [TU] we plot [ a/Fe] vs. time for our Case B model and the observed values from five 
red giants in Sextans (IShetrone et al.ll200ll ; iTolstoy et al.l 120031 1 . 



The evolution of [a/Fe] is very similar to that of [O/Fe]. Our models are generally 
consistent with the observations, but differ somewhat from the observations of the star S58. 
While our model predicts [a/Fe] rs —0.1 dex since 9 Gyr ago, the observed value for this 
one star is [a/Fe]~ —0.26. To fit the star S58, [a/Fe] should decrease abruptly as [Fe/H] 
increases. We tri ed to estimate th e age of S58 and the other four stars usi ng the new pho- 



tome try given by iLee et al.l (120031 ) and t he revised Yonsei-Y ale isochrones (jPemarque et al. 



2004 ) in the similar way to that used in ITolstoy et al.l (120031 ) . The estimated values for S58 
are 7 Gy r from (B — V ) color , and 6 Gyr from (V — /) color, which are similar to the value 
given by ITolstoy et al. The ages deriv ed for the other stars are larger than 14 Gyr, 

similar to those given by ITolstoy et al.l (120031 ) . 



In Figure[TT]we plot [a/Fe] versus [Fe/H] for the four radial regions in Sextans for Case B. 
[a/Fe] changes slowly about the value of 0.1, but then decreases sharply for [Fe/H]> —2. The 
predictions reach the solar value at [Fe/H]^ —1.7, and continue to drop to the the minimum 
value of ~ -0.1. The results from our modeling are consistent with the observations of four 
stars, but again differ for the star S58. 

Abundance ratios like [O/Fe] and [a/Fe] are affected by the rates of SN la and SN 
II during the evolution of a galaxy. In general SN II's have a significant effect from the 
beginning of the star formation history in galaxies with strong initial starbursts, while SN 
la's begin to play an important role in setting the abundance ratios after the time when the 
rate of SN la is a maximum, tsNTn (the typical tim e scale for the SN la maximum enrichment) 
(ITimmes et al.lll995l ; iMatteucci fc Recchi 1120011 ). The value of tsNia is as small as 0.3 — 0.5 
Gyr for elliptical galaxies and the bulges of spiral galaxies, in which the star formation rate 
was high in the beginning, while it is larger, 1.5 Gyr, for our Galaxy, and can be as large as 
7 — 8 Gyr for irregular galaxies, in which stars formed over a much more extended period of 
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time (IMatteucci h Recchi II2001I ). In Figure [12] we show the evolution of the SN II and SN la 
rates for our Case B model. The SN II rate was the maximum in the beginning, decreasing 
thereafter until the end of star formation. The SN la rate increased slowly from zero since 
the beginning, reaching the maximum 13-11 Gyr ago, when the star formation was active 
in the inner regions, and decreasing thereafter. Thus, tgNia is 2 to 4 Gyr for the Sextans 
dSph. This value is much larger than that for elliptical galaxies, but is between those of our 
Galaxy and irregular galaxies. 

Other abundance patterns, [X/Fe] vs. [Fe/H] , can also be useful in constraining the 
star formation and the chemical enrichment history (ITolstoy et al.1 120031 ) . In Figures [TBI we 
display the abundance ratio of se veral elements for our Case B model. The observations 
presented by lShetrone et al.l (120011 ) are plotted for comparison. The model abundance ratios 
are generally in good agreement with the observational ones. It is noted that the model 
values for [V/Fe] and [Cr/Fe] are marginally larger and smaller than the observation values, 
respectively, but the differences are at the level of one sigma. 

In Figures dU we s how the abundance of pres ent stable isotopes normalized with respect 
to the solar abundance ([Anders fc Grevesse Ill989l ) for our Case B model. The ratios of most 
of the heavy elements are between 0.1 to 0.001. 



4. Discussion 
4.1. Comparison to other models 



Ikuta fc Arimoto I (120021 ) investigated the star formation and chemical enrichment his- 
tory of three dSphs (Draco, Sextans, and Ursa Minor) using a chemical evolution model and 
a simulation program for CMDs. Like us, they adopted a Salpeter IMF and a closed box 
model. They concluded that the observed abundance of the giants and the CMD of the stars 
in Draco could be explained by very low star formation rates that continued for a long time 
since the beginning (3.9-6.5 Gyr). When they adopted much steeper IMFs (x = —1.75 to 
-2.15) and used a star formation rate similar to that of the solar neighborhood, they found 
that the star formation time could be shorter, 1.6-2.2. Gyr. They did not model Sextans, 
but noted that Sextan's SFH might be similar to that of Draco because of the analogous 
features in the CMD, a narrow RGB and red HB. Our result, that star formation continued 
for 8 Gyr in Sextans, is consistent with their inference based on their modeling of Draco 
with the same IMF. 



In contrast to our results, the star formation rate in the llkuta &: Arimoto I (120021 ) models 
does not decrease rapidly because it depends on the gas density, which declines slowly due 
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to the very low star formation rates. Hence, the resulting SFH should differ qualitatively 



from o urs. Such differences will be reflected in the evolution of [Fe/H]. In llkuta &: Arimoto 
( 120021 ) (see B and C models in their Figure 8), [Fe/H] increases slowly from the beginning 



reaching [Fe/H] = -2.4 to -2.7 dex at 1 Gyr, and then increases faster reaching [Fe/H] = -1.4, 
a few Gyr later. In our model, metallicity increases rapidly during the first 1 Gyr, but then 
changes little (Figures [7]). This behavior is c onsistent with the re sults for the first galaxies 
(similar to the Local Group dSphs) gi ven bvlKawata et all (120061) . The differences between 
the results presented here and those of Ikuta fc Arimoto ( 20021 ) are due to difference in the 
star formation rate. 



Lanfranchi fc Matteucci I (120041 ) studied the chemical evolution of Sextans and other 



dwarf spheroidal galaxies using the outflow model of chemical evolution, and found that 
their best model to match the observed abundance ratios of these dSphs is a combination 
of a high wind efficiency and a very low star formation efficiency. Their modeling results 
for [Si/Fe], [Mg/Fe], and [Ca/Fe] in the case of Sextans are consistent with the observed 
values (in their Figure 4), like ours. It is interesting that both our models based on a 
closed-box model and theirs based on a high galactic wind model are consistent with the 
observed abundance ratio. It is noted that we derived the star formation rate directly from 
the color-magnitude diagram, while they inferred it from the color-magnitude diagrams. 



4.2. Radial Gradients of Stellar Populations in dSphs 



Younger stellar populations with higher metallicity are often concentrated toward the 
centers of dSphs, a phenomenon often referred to as radial population gradients. For example, 
red h orizontal branch stars are more cen t rally concentrated th an blue horizontal branch 
stars (IHarbeck et al.ll200ll ; iLee et al.ll2003l ). lHarbeck et al.l (1200 ll ) discussed several possible 
explanations of these population gradients and concluded that metallicity is likely to be a 
main driver in some galaxies, inc luding Sextans, while age is the critical factor in others, 
such as Carina. iRizzi et al.l (120041 ) investigated the spatial distribution of the SFH of Carina, 
Sculptor, and Sextans and concluded that the driver was age for Carina, metallicity for 
Sculptor, and both of age and metallicity for Sextans. All of these require some spread in 
the characte ristics of the stel l ar pop ulations and both age and metallicity have been invoked 
for Sextans. iBellazzini et al.l (120011 ) argued that Sextans has two old populations of different 
metallicity, [Fe/ H] =-2 .3 and -1.8, but there is no other observational evidence for such a 



population mix (ILee et al.l 120031 ). Lastly, ILee et al.l (120031 ) concluded that the wide color 
distribution for the upper red giant branch is caused by continuous star formation, and that 
the metallicity increases during the period of a few Gyr, leading to A[Fe/H]?a 0.2 dex. 
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Our results support the findings of iLee et al.l (120031 ) and lRizzi et al.l (12004 ). The spatial 
variation of the star formation rate is obvious in this study: the star formation rate is higher 
in the central region than in the outer region. As star formation proceeds, the metallicity 
increases with time. This results in the radial gradient of metallicity, which is consistent 



with the radial gradient of the stellar population in Sextans ( jHarbeck et al.ll200ll ; ILee et al. 



20031 ). We conclude that the primary driver for the population gradient in Sextans is the 



SFH. It does not make any sense to argue for only either age or metallicity being the cause 
since they are related. 



4.3. Blue Straggler Stars 



If blue straggler stars are intermediate age main sequence stars, then they formed be- 
tween ~ 4 and ~ 7 Gyr ago. iLee et al.l (120031 ) pointed out, that in this case one expects 
red giant clump stars to be present between the red and blue horizontal branches, which are 
not seen in the observed CMD. Therefore, it is unlikely that these blu e strag glers are indeed 
intermediate-age normal main sequence stars. Similarly, ICarrea et al.l (120021 ) concluded that 
the U rsa Minor blue s tragg ler populations originates in the old population. On the other 



hand, iMapelli et al. I (120091 ) studied the spatial distribution of blue straggler stars based 
on wide field imaging, concluding that most of the blue straggler stars in the Fornax dSph 
are young main sequence stars, while the blue straggler stars in the Sculptor dSph are old 
binaries with mass-transfer. 

Our model cannot clarify the nature of the blue straggler stars, but can illustrate dif- 
ferences arising in the evolution of Sextans in the two scenarios. As we have shown, the 
SFH from 4 to 7 Gyr ago is where changes occur. However, the effect of this change on the 
chemical evolution is negligible because the number of the blue straggler stars is small in 
comparison to the rest of the population and the bulk of stars formed prior to this time. 



5. Summary and Conclusion 



Using the SMART model developed by Yuk & Lee (2007), we derive the SFH and 
chemical evolutio n of the Sextans dSph from the VI photometry of stars over the 42' x 28' 
field described by ILee et al.l (120031 ). Our principals results are summarized below. 



1. Independent of our assumption regarding the nature of blue straggler stars, the princi- 
pal star formation episode continued for several Gyr. From models in which we exclude 
the blue stragglers, we find that about 65 % of the stars formed within the first two 
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Gyr, about 25 % of stars formed between 13 Gyr ago and 11 Gyr ago, and the bulk 
of the remainder formed up until 7 Gyr ago. The models in which we include the blue 
stragglers and treat them as intermediate age main sequence stars, differ principally 
in that star formation continues at a low level until 4 Gyr ago. While the effect of the 
blue straggler stars on the SFH between 7 and 4 Gyr ago is noticeable, their effect on 
the chemical evolution is negligible. 

The mean metallicity of the stars increases rapidly up to [Fe/H]=-1.6 in the central 
region and [Fe/H]=-1.8 in the outer region for stars created within the first Gyr, and 
varies slowly for younger stars. 

Our derived abundance ratios of several elements are ge nerally in good agre ement 
with measurements based on high resolution spectroscopy (jShetrone et al.ll200ll ). Our 
derived abundance ratios for one star in the spectroscopic sample, S58, are much higher 
than the observed values. The reason for this difference is not known. 



There are significant differences in the SFH as a function of radius. Star formation was 
more vigorous and longer lived in the central region than in the outer region, resulting 
in a metallicity gradient. We conclude that the primary driver for the radial stellar 
population gradient in this galaxy is the SFH. 



Recently there are several works showing the importance of outflow (via galactic wind 
or via external mechanisms like ram pressure stripping) in the chemical e volution of dSphs 



flLanfranchi fc Matteucci 112004 120071 : iRecchi et al 



20081 : iMatteucci 1120081 ) . We are planning 



to investigate the effect of outflow in the chemical evolution of dSphs using the SMART model 
in the subsequent paper. 
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Table 1. Properties of the regions in the Sextans dSph 



Region a 


Radius b 


Area c 


iV 


iV 2 e 


Pi f 


p2 g 


Rl 


R < 7.'3 


167.4 


4318 


4259 


25.8 


25.4 


R2 


7.'3 < R < 11' 


212.7 


4697 


4633 


22.1 


21.8 


R3 


11' < R< 15' 


312.3 


5431 


5363 


17.4 


17.2 


R4 


15' < R 


483.6 


5864 


5795 


12.1 


12.0 



a The name of each region. 

b The boundary of each region in arcmin. 

c The area of each region in arcmin 2 . 

d Total number of stars in each region. 

e The number of stars excluding the blue straggler stars in 
the pentagon shown in Figure [TJ 

f Number density of stars in each region (= Ni / Area). 

g Number density of stars without blue stragglers in each 
region (= N 2 / Area). 



Table 2. Stars in the Sextans dSph with spectroscopic data by lShetrone et al.l (120011 ) 



Star ID a 


R.A. b 


Decl. c 


V d 


B — V c 


y-/ f 


Age g 


[Fe/H] h 


[a/Fe] 1 


S35 


10 h 13 m 57. s 72 


-1°38'53".0 


17.30 


1.45 


1.37 


15 


-1.93 ±0.11 


0.13 ±0.06 


S56 


10 h 12 m 41. s 86 


-1°45'27".0 


17.33 


1.43 


1.39 


15 


-1.93 ±0.11 


0.12 ±0.07 


S49 


10 h 13 m ll. s 57 


-1°43'01".4 


17.55 


1.16 


1.24 


15 


-2.85 ±0.13 


-0.01 ±0.11 


S58 


10 h 12 m 16. s 13 


-1°45'47".9 


17.65 


1.27 


1.29 


6 


-1.45 ±0.12 


-0.26 ±0.07 


S36 


10 h 13 m 04?96 


-1°39'13".4 


17.94 


1.16 


1.18 


15 


-2.19 ±0.12 


0.10 ±0.07 



a ID given bv lSuntzeff et alj (|l993h . Corresponding IDs in lLee et all |2003j) are 162, 169, 195, 212, and 255. 

b R. A. (2000) given bv lLee et af] [|2003h . 

c Dec. (2000) given bv lLee et all d2003h . 

d V magnitude given bv lLee et al. 

C B — V color given bv lLee et al. | l l2003h . 

{ V - I color given bv lLee et al] (|2003| V 

g Age of star in Gyr estimated by iTolstov et al. (2003) using the Yale-Yonsei isochrones jYi et al.ll200lh . New 
estimate for S58 is 7 Gyr based on the revised Yale-Yonsei isochrones dPemarque et al.ll2004F 

h [Fe/H] given bv lTolstov et alj j2003h . 

'[a/Fe] given bv lTolstov et alj j2003h . 
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Table 3. Model parameters 



item 


value 


distance modulus (m — M) a 


19.67 


reddening E(B — V) b 


0.03 


IMF index c 


-1.35 


mass range 


O.O8Af - 40M o 


chemical evolution model 


closed box model 


initial composition d 


primordial value 



a Determined in this study. 
b Determined in this study. 



c In dex of the Salpeter's IMF (x) (ISalpeter 
1953). 



d Homogeneous big 
JWalker et alJll99lh . 



bang composition 
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Table 4. Star formation history for Case A 



Age (Gyr) a 


SFR (Rl) b 


SFR (R2) c 


SFR (R3) d 


SFR (R4) c 





- 4 


0.00 ±0.01 


0.00 ±0.01 


0.00 ±0.01 


0.00 ±0.01 


4 


- 7 


0.03 ±0.01 


0.02 ±0.01 


0.01 ±0.01 


0.01 ±0.01 


7 


- 9 


0.10 ±0.05 


0.07 ±0.04 


0.06 ±0.03 


0.04 ±0.01 


9 


- 11 


0.19 ±0.17 


0.16 ±0.12 


0.14 ±0.08 


0.07 ±0.04 


11 


- 13 


0.61 ±0.26 


0.53 ±0.19 


0.30 ±0.13 


0.13 ±0.07 


13 


- 14 


0.99 ±0.34 


0.86 ±0.25 


0.78 ±0.18 


0.64 ±0.11 


14 


- 15 


0.99 ±0.14 


0.87 ±0.12 


0.78 ±0.10 


0.66 ±0.06 


a Ag 


;e bin in 


Gyr. 









b ' c ' d ' e The star formation rates of the central region Rl (R < 7/3), 
the inner region R2 (7/3 < R < 11'), the intermediate region R3 
(11' < R < 15'), and the outer region R4 (15' < R), respectively. 
These were normalized with respect to the lifetime average star 
formation rate for the entire observed region of Sextans, 1.7 x 10~ 4 
M yr" 1 . 
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Table 5. Star formation history for Case B 



Age (Gyr) a 


SFR (Rl) b 


SFR (R2) c 


SFR (R3) d 


SFR (R4) c 





- 4 


0.00 ±0.01 


0.00 ±0.01 


0.00 ±0.01 


0.00 ±0.01 


4 


- 7 


0.00 ±0.01 


0.00 ±0.01 


0.00 ±0.01 


0.00 ±0.01 


7 


- 9 


0.05 ±0.02 


0.04 ±0.01 


0.02 ±0.01 


0.02 ±0.01 


9 


- 11 


0.13 ±0.09 


0.12 ±0.06 


0.08 ±0.04 


0.05 ±0.02 


11 


- 13 


0.65 ±0.19 


0.41 ±0.14 


0.29 ±0.09 


0.10 ±0.04 


13 


- 14 


1.07 ±0.29 


1.06 ±0.22 


0.88 ±0.15 


0.68 ±0.09 


14 


- 15 


1.07 ± 0.11 


1.06 ±0.10 


0.89 ±0.08 


0.77 ±0.06 


a Ag 


;e bin in 


Gyr. 









b ' c ' d The star formation rates of the central region Rl (R < 7.'3), 
the inner region R2 (7.'3 < R < 11'), the intermediate region R3 
(11' < R < 15'), and the outer region R4 (15' < R), respectively. 
These were normalized with respect to the lifetime average star 
formation rate for the entire observed region of Sextans, 1.7 x 10 -4 
M yr" 1 . 
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Table 6. Mean metallicity of each region in the Sextans dSph 



Region 


< Fe/H > (Case A) 


< Fc/H > (Case B) 


Rl 


-1.69 


-1.68 


R2 


-1.76 


-1.73 


R3 


-1.86 


-1.83 


R4 


-1.99 


-1.97 


Entire field 


-1.78 


-1.76 
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Fig. 1. — I — (V — I) diagrams of the four regions at different radii from the center of the 
Se xtans dSph. Op e n circ les represent the stars with high resolution spectroscopic data given 
by lShetrone et al.l (1200 ll ). Pentagons represent the region of the blue straggler stars. 
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0.5 1 1.5 

(V-I) 




0.5 1 1.5 



(v-i) 

Fig. 2. — The observational CMD (upper panel) and the best fit synthetic CMD for Case A 
(lower panel) of the Sextans dSph. The solid line corresponds to the 60% completeness level 
of the photometry. 
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Fig. 3. — Star formation history of the four regions in the Sextans dSph for Case A where blue 
straggler stars were considered as main-sequence stars with intermediate age for modeling. 
The star formation rates were normalized with respect to the lifetime average star formation 
rate for the entire observed region of Sextans, 1.7 x M Q yr _1 . The vertical bars indicate 
the asymmetric errors of 1 a level. Note that a primary star formation continued for long 
until 7 Gyr ago, followed by a minor star formation from 7 to 4 Gyr ago 
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Fig. 4. — Radial distribution of the star formation rates of the Sextans dSph for Case A. 
The star formation rates were normalized with respect to the lifetime average star formation 
rate for the entire observed region of Sextans, 1.7 x 10~ 4 M & yr _1 . The vertical bars indicate 
the symmetric errors of 1 a level. Note the radial gradient of the star formation rate for 
extended period of time. 
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Fig. 5. — Star formation history of the Sextans dSph for Case B where blue straggler stars 
were not considered for modeling. The star formation rates were normalized with respect to 
the lifetime average star formation rate for the entire observed region of Sextans, 1.7 x 1CT 4 
M & yr _1 . The vertical bars indicate the asymmetric errors of 1 a level. Note that a primary 
star formation continued for long until 7 Gyr ago. 
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Fig. 6. — Radial variation of the star formation rates of the Sextans dSph for Case B. The 
star formation rates were normalized with respect to the lifetime average star formation rate 
for the entire observed region of Sextans, 1.7 x 1CT 4 M & yr _1 . The vertical bars indicate 
the asymmetric errors of 1 a level. Note the radial gradient of the star formation rate for 
extended period of time. 
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0.01 0.1 1 10 

15 - Age [Gyr] 

Fig. 7. — Metal enrichment history for the four regions in the Sextans dSph for Case B. The 
higher the lines are, the closer the corresponding regions are to the galaxy center, (solid line 
for Rl, short-dashed line for R2, long-dashed line for R3, and dot-dashed line for R4). 
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Fig. 8. — Radial variation of the metallicity [Fe/H] of the Sextans dSph for Case B. The 
numbers beside lines represent the age of stars in Gyr. 
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Fig. 9. — Evolution of [O/Fe] for the four regions in the Sextans dSph Case B. 
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Fig. 10. — Evolution of [ot/Fe] for the four regions in the Sextans dSph for Case B derived in 
this study (solid lines), where [a/Fe] = ([Mg/Fe] + [ Ca/Fel + [Ti/FeD/3. Filled circles with 
error bars represent the observational data given by IShetrone et al.l (120011 ) 
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Fig. 11. — [a/Fe] ratios versus [Fe/H] for the four regions in the Sext ans dSph for Case B . 
Filled circles with error bars represent the observational data given by lTolstoy et al.l (120031 ). 
and the lines represent the values derived in this study. 
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Fig. 12. — Evolution of the SN la (lower group) and SN II (upper group) rate for the four 
regions in the Sextans dSph for Case B. 
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Fig. 13. — Abundance ratios for several elements of the Sexta ns clSph for Ca s e B. Filled 
circles with error bars represent the observational data given by lShetrone et al.l (120011 ). and 
the solid lines represent the results derived in this study. 
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Fig. 14. — The present stable isotopes of the Sextans dSph normalized to the solar abundance 
for Case B. The starlets represent the isotopes most abundant in the Sun. 



